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We discuss the possible sources of dark matter axions in the early universe. In the standard thermal scenario, 
an axion string network forms at the Peccei-Quinn phase transition T ~ / a and then radiatively decays into a 
cosmological background of axions; to be the dark matter, these axions must have a mass m a ~ 100 [ieV with 
specified large uncertainties. An inflationary phase with a reheat temperature below the PQ-scale Tleh 1$ /a can 
also produce axion strings through quantum fluctuations, provided that the Hubble parameter during inflation is 
large Hi<l/ a ; this case again implies a dark matter axion mass m a ~ 100 peV. For a smaller Hubble parameter 
during inflation HiiS/a, ‘anthropic tuning’ allows dark matter axions to have any mass in a huge range below 
m a < ImeV. 


1. INTRODUCTION 


2. STANDARD AXION COSMOLOGY 


The axion has remained a popular dark matter 
candidate because of its enduring motivation as 
an elegant solution to the strong CP-problenr [jj]. 
Despite early hopes of discovery, it turned out 
that in order to be consistent with accelerator 
searches and astrophysics, the axion must be 
nearly ‘invisible’ and extremely light. Its cou¬ 
plings and mass are inversely proportional to the 
(large) Peccei-Quinn scale / a as in 


m a = 6.2 x 10 _6 eV 


^ 10 12 GeV ^ 


( 1 ) 


Accelerator constraints have been largely super¬ 
seded by those from astrophysics; because the 
axion is so weakly coupled, volume effects can 
compete with other surface and convective stel¬ 
lar energy loss mechanisms. The strongest as- 
trophysical constraints on the axion mass de¬ 
rive from studies of supernova 1987a and con¬ 
servative estimates yield m a ^ 10meV||. The 
present programme of large-scale axion search 
experiments [|| are sensitive to a mass range 
m a ~ 1-10 /ieV, which has been chosen for a va¬ 
riety of historical and technological reasons. Our 
primary focus here, however, is not on constraints 
on the viable axion mass range, but rather on ef¬ 
forts to predict the mass of a dark matter axion 
from cosmology. 


The cosmology of the axion is determined by the 
two energy scales / a and Aqcd- The first im¬ 
portant event is the Peccei-Quinn phase tran¬ 
sition which is broken at a high temperature 
T ~ / a > 10 9 GeV. This creates the axion, at this 
stage an effectively massless pseudo- Goldstone 
boson, as well as a network of axion strings j| 
which decays gradually into a background cosmic 
axions ||. (Note that one can engineer models in 
which an inflationary epoch interferes with the ef¬ 
fects of the Peccei-Quinn phase transition, as we 
shall discuss in the next section.) At a much lower 
temperature T ~ Aqcd after axion and string for¬ 
mation, instanton effects ‘switch on’, the axions 
acquire a small mass, domain walls form || be¬ 
tween the strings JJ] and the complex hybrid net¬ 
work annihilates into axions in about one Hubble 
time JJ. 

There are three possible mechanisms by which 
axions are produced in the ‘standard thermal sce¬ 
nario’: (i) thermal production, (ii) axion string 
radiation and (iii) hybrid defect annihilation 
when T = A^qcd- Axions consistent with the as- 
trophysical bounds must decouple from thermal 
equilibrium very early; their subsequent history 
and number density is analogous to the decoupled 
neutrino, except that unlike a lOOeV massive neu¬ 
trino, thermal axions cannot hope to dominate 
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the universe with m a lOmeV. We now turn to 
the two dominant axion production mechanisms, 
but first we address an important historical di¬ 
gression. 

2.1. Misalignment misconceptions 

The original papers on axions suggested that 
axion production primarily occurred, not through 
the above mechanisms, but instead by ‘misalign¬ 
ment’ effects at the QCD phase transition ||]. 
Before the axion mass ‘switches on’, the axion 
field 9 takes random values throughout space in 
the range 0 to 27 t; it is the phase of the PQ- 
field lying at the bottom of a U{ 1) ‘Mexican 
hat’ potential. However, afterwards the potential 
becomes tilted and the true minimum becomes 
9 = 0, so the field in the ‘misalignment’ pic¬ 
ture begins to coherently oscillate about this min¬ 
imum; this homogeneous mode corresponds to the 
‘creation’ of zero momentum axions. Given an 
initial rms value 9 X for these oscillations, it is rela¬ 
tively straightforward to estimate the total energy 
density in zero momentum axions and compare 
these to the present mass density of the universe 
(assuming a flat fl = 1 FRW cosmology) 

n a , hom » 2A/r- 2 0 I 2 f(0 i )f^-^'j (2) 

where A « 3 ±1 accounts for both model- 
dependent axion uncertainties and those due to 
the nature of the QCD phase transition, and h 
is the rescaled Hubble parameter. The function 
f(0) is an anharmonic correction for fields near the 
top of the potential close to unstable equilibrium 
9 ss 7r, that is, with f(0) = 1 at the base 9 « 0 
and diverging logarithmically for 9 —► 7r jn|. If 
valid, the estimate (J2j) would imply a constraint 
m a > 5 geV for the anticipated thermal initial con¬ 
ditions with 9\ = 0(1) |,|. 

As applied to the thermal scenario, the expres¬ 
sion @) is actually a very considerable underes¬ 
timate for at least two reasons: First, the ax¬ 
ions are not ‘created’ by the mass ‘switch on’ at 
t = tQCDi they are already there with a specific 
momentum spectrum g(k) determined by dynam¬ 
ical mechanisms prior to this time. The actual 
axion number obtained from g{k) is much larger 
than the rms average assumed in (|j) which ig¬ 


nores the true particle content. Secondly, this 
estimate was derived before much stronger topo¬ 
logical effects were realized, notably the presence 
of axion strings and domain walls. In any case, 
these nonlinear effects complicate the oscillatory 
behaviour considerably, implying that the homo¬ 
geneous estimate (j{2|) is poorly motivated. 

2.2. Axion string network decay 

Axions and axion strings are inextricably inter¬ 
twined. Like ordinary superconductors or super¬ 
fluid 4 He, axion models contain a broken U( 1)- 
symmetry and so there exist vortex-line solutions. 
Combine this fact with the Peccei-Quinn phase 
transition, which means the field is uncorrelated 
beyond the horizon, and a random network of 
axion strings must inevitably form. An axion 
string corresponds to a non-trivial winding from 
0 to 27 t of the axion field 9 around the bottom 
of its ‘Mexican hat’ potential. It is a global 
string with long-range fields, so its energy per 
unit length g has a logarithmic divergence which 
is cut-off by the string curvature radius RtSt, that 
is, g ss 2wf£ ln(t/6 ), where the string core width 
is S « f- 1 . The axion string, despite this loga¬ 
rithmic divergence, is a strongly localized object; 
if we have a string stretching across the horizon 
at the QCD temperature, then In {t/5) ~ 65 and 
over 95% of its energy lies within a tight cylinder 
enclosing only 0.1% of the horizon volume. To 
first order, then, the string behaves like a local 
cosmic string, a fact that can be established by 
a precise analytic derivation and careful compar¬ 
ison with numerical simulations [ |lT[ . 

After formation and a short period of damped 
evolution, the axion string network will evolve to¬ 
wards a scale-invariant regime with a fixed num¬ 
ber of strings crossing each horizon volume (for 
a cosmic string review see ref. @)- This grad¬ 
ual demise of the network is achieved by the pro¬ 
duction of small loops which oscillate relativis- 
tically and radiate primarily into axions. The 
overall density of strings splits neatly into two 
distinct parts, long strings with length £ > t 
and a population of small loops £ < t, that is, 
P = Poo + Pl- High resolution numerical simula¬ 
tions confirm this picture of string evolution and 
suggest that the long string density during the ra- 
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diation era is ~ 13 p,/t 1 2 Q. To date, analytic 
descriptions of the loop distribution have used 
the well-known string ‘one scale’ model, which 
predicts a number density of loops defined as 
p£n(£,t) d£ = pL{£,t)d£ in the interval £ to £ + d£ 
to be given by 


4a 1//2 (l + «;/a) 3 / 2 
(£ + nt) 5 / 2 t 3 / 2 


(3) 


where a is the typical loop creation size relative 
to the horizon and k ss 65/[27rln {t/5)\ is the loop 
radiation rate |Q . Once formed at t = to with 
length a typical loop shrinks linearly as it de¬ 
cays into axions £ = £q — n{t — to)- The key un¬ 
certainty in this treatment is the loop creation 
size a, but compelling heuristic arguments place 
it near the radiative backreaction scale, a ~ k. 
(If this is the case, we note that the loop contri¬ 
bution is over an order of magnitude larger than 
direct axion radiation from long strings.) 

String loops oscillate with a period T = £/2 and 
radiate into harmonics of this frequency (labelled 
by n), just like other classical sources. Unless a 
loop has a particularly degenerate trajectory, it 
will have a radiation spectrum P„ oc n~ q with a 
spectral index q > 4/3, that is, the spectrum is 
dominated by the lowest available modes.Q Given 
the loop density (||), we can then calculate the 
spectral number density of axions dn a /<iw, which 
turns out to be essentially independent of the ex¬ 
act loop radiation spectrum for q > 4/3. From 
this expression we can integrate over to to find 
the total axion number at the time £qcd, that 
is, when the axion mass ‘switches on’ and the 
string network annihilates. Subsequently, the ax¬ 
ion number will be conserved, so we can find the 
number-to-entropy ratio and project forward to 
the present day. Multiplying the present number 
density by the axion mass m a yields the over¬ 
all axion string contribution to the density of the 
universe & 

/ lO _6 eV\ 118 

Ua,string ~ 110A h~ 2 [ - /(«/«) , (4) 

V m a / 

1 Historically, there has been some debate on the radiation 

spectrum issue but, the reader is referred elsewhere for 

further details ||lq|. 


where 

/(«/«) 



(5) 


The key additional uncertainty from the string 
model is the ratio a/t t ~ 0(1), which should 
be clearly distinguished from particle physics and 
cosmological uncertainties inherent in A and h 
(which appear in all estimates of fl a ). With a 
Hubble parameter near h = 0.5, the string esti¬ 
mate ([I|) tends to favour a dark matter axion with 
a mass m a ~ lOO^eV, as we shall discuss in the 
conclusion. A comparison with ([/[) confirms that 
string is well over an order of magnitude larger 
than the ‘misalignment’ contribution. 


2.3. Hybrid defect annihilation 

Near the QCD phase transition the axion ac¬ 
quires a mass and network evolution alters dra¬ 
matically because domain walls form. Large field 
variations around the strings collapse into these 
domain walls, which subsequently begin to dom¬ 
inate over the string dynamics. This occurs 
when the wall surface tension a becomes com¬ 
parable to the string tension due to the typical 
curvature a ~ p/t. The demise of the hybrid 
string-wall network proceeds rapidly, as demon¬ 
strated numerically |/J. The strings frequently in¬ 
tersect and intercommute with the walls, effec¬ 
tively ‘slicing up’ the network into small oscillat¬ 
ing walls bounded by string loops. Multiple self¬ 
intersections will reduce these pieces in size until 
the strings dominate the dynamics again and de¬ 
cay continues through axion emission. 

An order-of-magnitude estimate of the demise 
of the string-domain wall network indicates that 
there is an additional contribution ||l6|| 

/in-6 P V\ 118 

H a , dw ~ 0(10) Ah~ 2 - . (6) 

\ m a / 


This ‘domain wall’ contribution is ultimately due 
to loops which are created at the time ~ tQCD- 
Although the resulting loop density will be similar 
to (|^), there is not the same accumulation from 
early times, so it is likely to be subdominant (id) 
relative to !)• More recent work, (l7| questions 
this picture by suggesting that the walls stretch¬ 
ing between long strings dominate and will pro- 
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duce a contribution anywhere in the wide range 
^a,dw ^ (1 44)r^ ajS tring? however, this assertion 
requires stronger quantitative support. Overall, 
like most effects]^] the domain wall contribution 
will serve to further strengthen the string bound 
([|) on the axion. 

Up to this point we have only considered the 
simplest axion models with a unique vacuum 
N = 1, so what happens when N > 1? In this 
case, any strings present become attached to N 
domain walls at the QCD-scale. Such a network 
‘scales’ rather than annihilates, and so it is cos- 
mologically disastrous being incompatible (at the 
very least) with CMB isotropy. 

3. INFLATIONARY AXION MODELS 

The relationship between inflation and dark mat¬ 
ter axions is rather mysterious. Its significance 
depends on the magnitude of the Peccei-Quinn 
scale / a relative to two key inflationary param¬ 
eters, (i) the reheat temperature of the universe 
T re h at the end of inflation and (ii) the Hubble 
parameter Hi as the observed universe first exits 
the horizon during inflation. Inflation is irrele¬ 
vant to the axion if T re h ^ / a because, in this case, 
the PQ-symmetry is restored and the universe re¬ 
turns to the ‘standard thermal scenario’ in which 
axion strings form and the estimate (|j) pertains. 
Consider, then, the two inflationary axion scenar¬ 
ios with T reh <f a . 

3.1. Anthropic misalignment and quantum 
fluctuations (Case 1: Hi < fj 

In an inflationary model for which / a > Hi > 
T re h, the 0-parameter or axion angle will be set 
homogeneously over large inflationary domains 
before inflation finishes |l8|]. In this case, the 
whole observable universe emerges from a single 
Hubble volume in which this parameter has some 
fixed initial value 0p Because the axion remains 
out of thermal equilibrium for large / a , subse¬ 

2 We note briefly that it is also possible to weaken any 
axion mass bound through catastrophic entropy produc¬ 
tion between the QCD-scale and nucleosynthesis, that is, 
in the timescale range 10“ 4 si$t e nt 10 —2 s. Usually this 
involves the energy density of the universe becoming tem¬ 
porarily dominated by an exotic massive particle with a 
tuned decay timescale. 


quent evolution and reheating does not disturb 8, 
until the axion mass ‘switches on’ at T ~ Aqcd- 
Afterwards, the field begins to oscillate coher¬ 
ently, because it is misaligned by the angle 0i 
from the true minimum 0 = 0. This homoge¬ 
neous mode corresponds to a background of zero 
momentum axions and it is the one circumstance 
under which the misalignment formula (|2j) actu¬ 
ally gives an accurate estimate of the relative ax¬ 
ion density fl a . 

By considering the dependence ff a h oc 6 ( in 
(f|), we see that inflation models have an intrinsic 
arbitrariness given by the different random mag¬ 
nitudes of 0; in different inflationary domains | jl8| . 
While a large value of / a » 10 12 GeV might have 
been thought to be observationally excluded, it 
can actually be accommodated in domains where 
0i « 1. This may seem highly unlikely but, if 
we consider an infinite inflationary manifold or 
a multiple universe scenario including ‘all possi¬ 
ble worlds’, then life as we know it would be ex¬ 
cluded from those domains with large 0, = 0(1) 
because the baryon-to-axion ratio would be too 
low 1191. Thus, accepting this anthropic selection 
effect, we have to concede that axions could be 
the dark matter f2 a ps 1 if we live in a domain 
with a ‘tuned’ 0-parameter|: 


0i 


0.3 A-V 2 3 /* 

110- 6 eV/ 


0.6 


(7) 


For 0j ps 0(1), this suggests an axion with m a ~ 
5peV (h = 0.5), though actually inflation makes 
no definite prediction from (|fj) beyond apparently 
specifying m a > 10 peV. But even this restriction 
is not valid; if we observe (2!) carefully we see 
that we can also obtain a dark matter axion for 
higher m a by fine-tuning 0; near 7r. The anhar- 
monic term f(0) with an apparent logarithmic di¬ 
vergence allows S2 a ss 1 for a much heavier dark 
matter axion ©■ 

This simple ‘anthropic tuning’ picture is sig¬ 
nificantly altered by quantum effects. Like any 
minimally coupled massless field during inflation, 
the axion will have a spectrum of quantum ex¬ 
citations associated with the Gibbons-Hawking 
temperature T ~ H/2n. This implies the field 

3 This is not quite in the spirit of the original motivation 

for the axion! 
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will acquire fluctuations about its mean value 0; 
of magnitude 88 = H/27r/ a , giving an effective 
rms value 0 2 ff = (0i + 89) 2 . Even if our uni¬ 
verse began in an inflationary domain with 0i = 0, 
there will be a minimum misalignment angle set 
by 88; this implies that we cannot always fine- 
tune 8i in 0) such that f2 Q <l. Worse still, the 
fluctuations 88 imply isocurvature fluctuations 
in the axion density Q given by (|j), that is, 
8 oc 0 2 ff « 8f + 29i88 + 89' 2 . Such density fluc¬ 
tuations will also create cosmic microwave back¬ 
ground (CMB) anisotropies which are strongly 
constrained, ST/T ~ 8p/p ~ 288/9 e g< 10 -5 . 
Combining the quantum fluctuation 86 with the 
0-requirement for a dark matter axion ( 0 ), we ob¬ 
tain a strong constraint on the Hubble parameter 
during inflation {§ 


Hi < f0 9 GeV 


10~ 6 eV \ °' 4 

m a J 


( 8 ) 


which is valid for small masses m a 10 peV. 
Here, Hi is the Hubble parameter as fluctuations 
associated with the time when the microwave 
anisotropies first leave the horizon some 50-60 e- 
foldings before the end of inflation. We conclude 
from (0) that inflation in the / a > Hi > T re h 
regime must have a small Hubble parameter Hi 
or the dark matter axion must be extremely light, 
m a <K 1 peV. Inflation now has many guises and 
models exist with Hubble parameters anywhere 
in the range Hi ~ 10 2 10 14 GeV, so a detectable 
inflationary axion with m a ~ 1 peV is possible in 
principle. Note, however, that the constraint (0) 
can be circumvented in more complicated multi- 
field inflation models p^| . 

Quantum effects also constrain the ‘anthropic 
tuning’ 0; —> 7T which is required for a larger 
mass axion m a >10^eV. For simple inflation 
models, excessive isocurvature fluctuations imply 
the dark matter axion mass is bounded above by 
m a <lmeV for Hi > 10 2 GeV |uj. Note that this 
is a considerably heavier dark matter axion than 
the oft-quoted inflationary limit m a < 5 peV. 


3.2. Axion string network creation during 
inflation (Case 2: H \ 

Even for a low reheat temperature T re h < / a , 
one can envisage inflation models with a large 


Hubble parameter during inflation Hi>/ a (such 
as chaotic inflation with Hi ~ 10 13 — 10 14 GeV). 
The quantum fluctuations in this case are suffi¬ 
cient to take the Peccei-Quinn field over the top 
of the potential leaving large spatial variations 
and topologically non-trivial windings [^|. We 
can interpret this as the Gibbons-Hawking tem¬ 
perature ‘restoring’ the PQ-symmetry Tcm^/a- 
As inflation draws to a close and H falls below 
/ a , these fluctuations will become negligible and 
a string network will form. Provided inflation 
does not continue beyond this point for more than 
about another 30 e-foldings, we will effectively re¬ 
turn to the ‘standard thermal’ scenario in which 
axions are produced by a decaying string net¬ 
work. So such low reheat inflation models are 
again only compatible with a dark matter axion, 
m a ~ 100/reV. 

We note that there is also a borderline scenario 
with Hi ~ fa and in which domain walls form but 
few strings. Since strings are required to remove 
them, such domian walls will be cosmologically 
unacceptable |m) . 

4. CONCLUSIONS 

We have endeavoured to provide an overview 
of axion cosmology focussing on the mass of a 
dark matter axion. First, the cosmological axion 
density was considered in the standard thermal 
scenario where the dominant contribution comes 
from axion strings. In this case there is, in prin¬ 
ciple, a well-defined calculational method to pre¬ 
cisely predict the mass m a of a dark matter ax¬ 
ion. For the currently favoured value of the Hub¬ 
ble parameter ( H 0 ss 60kms _1 Mpc _1 ), the esti¬ 
mate (0) predicts a dark matter axion of mass 
Tiia ~ 200 peV , where significant uncertainties 
from all sources approach an order of magnitude. 
The additional uncertainty in this string calcula¬ 
tion is the parameter ratio a/n ~ 1, that is, the 
ratio of the loop size to radiation backreaction 
scale. 

Secondly, we have reviewed inflationary axion 
cosmology showing that (i) many inflation mod¬ 
els return us to the standard thermal scenario 
with iria ~ 100/ieV, (ii) some inflation models are 
essentially incompatible with a detectable dark 
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matter axion and (iii), because of the possibility 
of ‘anthropic fine-tuning’, other inflation models 
can be constructed which incorporate a dark mat¬ 
ter axion mass anywhere below m a 1$ lrneV. 

We conclude that, while a dark matter ax¬ 
ion might possibly lurk anywhere in an enor¬ 
mous mass range below m a ^51 meV, the best- 
motivated mass for future axion searches lies near 
?7i a ~ 100/zeV, a standard thermal scenario pre¬ 
diction which is also compatible with a broad 
class of inflationary models. 
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